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Abstract. We revisit a well-studied solar flare whose X- 
ray emission originating from a simple loop structure was 
observed by most of the instruments on board SMM on 
November 12 1980. The X-ray emission of this flare, as ob- 
served with the XRP, was successfully modeled previously. 
Here we include a detailed modeling of the transition re- 
gion and we compare the hydrodynamic results with the 
UVSP observations in two EUV lines, measured in areas 
smaller than the XRP rasters, covering only some portions 
of the flaring loop (the top and the foot-points) . The single 
loop hydrodynamic model, which fits well the evolution of 
coronal lines (those observed with the XRP and the Fe xxi 
1354.1 A line observed with the UVSP) fails to model the 
flux level and evolution of the Ov 1371.3 A line. 

Key words: corona, transition region, flares, hydrody- 
namics 



1. Introduction 

Solar flares are very complex phenomena. They emit in 
a wide range of wavelengths, including radio, optical, UV 
and X-rays. They involve many physical effects such as, 
for instance, chromospheric evaporation, magnetic field re- 
connection and material ejection. 

There have been many examples of hydrodynamic 
models of flaring plasma confined in solar coronal loops: 
for a list of such models and some recent developments see 
Peres (2000")| . In this context [Peres et al. (1987)| modeled 
the X-ray emission of the well-studied solar fiarc, which 
occurred on November 12 1980 at 17:00UT. In particular, 
they used the Palermo-Harvard hydrodynamic loop model 
(Peres et al. (1982), thereafter PH), to give an interpreta- 
tion of the phenomena involved in this event. The mod- 
eled light curves were successfully compared with those 



observed in some X-ray lines with the XRP on board SMM 
in a raster region covering the flaring coronal loop. 

The numerical resolution in the transition region (TR), 
however, was not sufficient for a proper comparison with 
the EUV observations. Most EUV lines, in fact, are formed 
at temperat ures below lO^K. U sing the new version of 
the PH code (Betta et al. (1997)), having appropriate spa- 
tial resolution, here we revisit the modeling of this well- 
studied flare. In section 2 we describe the hydrodynamic 
loop model; in section 3 we introduce our new simulations; 
in section 4 we compare the PH code numerical calcula- 
tions with the observed line light curves; in section 5 we 
summarize our conclusions. 



2. The flaring loop model 
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The PH code ([Peres et al. (1982)D solves the one-fluid 
time-dependent differential equations of mass, momentum 
and energy conservation for the plasma confined in a solar 
magnetic loop, using a finite-difference numerical scheme. 
The model is one-dimensional, since the plasma is con- 
fined inside a semicircular flux tube where bulk motion 
and heat flow occur only along the magnetic field lines. 
The present version of the PH code ( Betta et al. (199"7)| ) 
can solve an asymmetric loop (e.g. pR.eale et al. (2000) ), 
however MacNeice et al. (1985) inferred that the fiare X- 
ray emission of the Nov 12 1980 fiare was symmetric re- 
spect to the apex, and therefore we limit to model half a 
loop. We also assume ionization and thermal equilibrium 
between ions and electrons and that the fiux tube keeps 
its geometric shape. The equations are solved with proper 
accuracy along the loop and during the entire fiare evolu- 
tion, since a regridding algorithm adapts the spatial grid 
whenever and wherever needed. The adaptive grid allows 
to reach very high spatial resolutions (as small as Im) 
using a variable number of grid points (~ 300 — 600). 

We consider the same loop parameters used in 
Peres et al. (1987), which were inferred from an accurate 



analysis of the Ha and X-ray images. 
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The flare event is simulated by switching on a very 
strong impulsive heating Q. We use the formulation 

Q{s,t)=Hs+Ho fit) his) 

where s is the coordinate parallel to the magnetic field 
lines, t is time, iJ^, uniform and constant, balances or- 
dinary coronal losses, while the second term in the right 
handside represents the impulsive heating. This term may 
be due to current dissipation. Its spatial and temporal de- 
pendencies are separated into two factors. In general we 
allow for various functional forms for fit) and his). Here 
we assume that the spatial term /i(s) has a Gaussian form 
of half- width a and center sq] fit) is a constant function 
for the first 180s and then decays exponentially with an 
e-folding time of 60s. 

We repeat the following simulations: 

A) heating applied at the top of the loop; 

B) heating applied in a small region near the foot-points. 

We compute the line emission from the temperature 
and density of the plasma along the loop obtained 
with the PH code, using the software ASAP written by 



Maggio et al. (1994) in IDL 



3. Numerical results 

The hydrodynamics of the coronal plasma subject to a 
strong flaring heating has already been described exten- 
sively in previous papers in the literature. Here we dwell 
on the most important details. 

The evolution of the temperature and the density at 
the top of the loop are shown in Fig. |l]. In the first 
minute of evolution, the coronal temperature increases 
from roughly 3 to 21 MK, for the model in which the 
energy is applied at the loop top (simulation A), and to 
around 15 MK for the model in which the energy is ap- 
plied at the base of the TR (simulation B). During the 
early flare we find the largest differences between the two 
models. 

In both simulations A and B, the coronal density in- 
creases while the impulsive heating is on. After 200s of 
the simulations, the hydrogen density at the top of the 
loop is of the order of lO^^m"^. In the pre-flare conditions 
the plasma hydrogen density at the top of the loop was 
instead ~ 7 lO^^m"^. 

As the coronal temperature increases, the conductive 
front propagates through the TR (the latter at the same 
time becomes steeper and steeper), and the chromospheric 
plasma warms up and, after ~ 15s, expands rapidly toward 
the corona ("chromospheric evaporation"). 

In simulation A, velocities at temperatures above lO^K 
are always toward the loop top after the first 24s and until 
the impulsive heating is active, i.e until 360s. In ^ 30s 
since the beginning of this simulation, the plasma velocity 
is ^ 300km s~^ in the corona. On the other hand, when 
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Fig. 1. Temporal evolution of the temperature and hydro- 
gen density at the top of the loop for the two models. 

the heating is applied at the base of the loop (simulation 
B), the plasma is pushed up since the very beginning and 
reaches velocities ~ 400km s~^ after 20s. 

In both simulations a descendant flow of material oc- 
curs in the chromosphere, below the region being ablated 
and evaporated; there the plasma pressure is much lower 
than the pressure of the plasma evaporated in the corona 
and the material is shocked and compressed toward the 
bottom; velocities have much lower values than coronal 
ones. This process is generally called " chromospheric con- 
densation" ( [Fisher 1987|) . 

As the plasma starts to evaporate, in both simulations 
the plasma temperature progressively increases along the 
whole loop; the emission measure at any temperature in- 
creases (as well as the EUV flux). The TR moves along 
the loop toward the foot-points during the whole heat- 
ing phase. After 240s, when the impulsive heating has de- 
creased to around 1/3 of its maximum value, the corona 
starts to cool and the TR moves gradually back to its 
initial position. 
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The simulations have been run to describe the plasma 
evolution for many minutes after the end of the heating. 
The loop cools and t he coronal density decreases until 
a thermal instability ( Field 1965 ) occurs: then the whole 
loop cools very rapidly and reaches typical chromospheric 
values (~ 10*K). As a uniform and constant heating term 
(Hs) is present in the equations, the whole loop returns 
approximately to its stationary pre-flare equilibrium after 
^ 3000s. The entire cycle - until the whole coronal at- 
mosphere returns exactly to the hydrostatic pre-flare con- 
ditions and the velocities decrease to negligible values - 
takes a few hours. 



4. Light curves 



This flare has been described by MacNeice et al. (1985) 
and pheng fc Pallavicini (1988)| (CP88, thereafter). We 
now summarize the most important features for a better 
comprehension of what is discussed in the following. 

It started at 17:00 UT and lasted less than 20 minutes. 
Four of the five instruments on board SMM registered 
the event: the Hard X-ray Burst Spectrometer (HXRBS), 
the Hard X-ray Imaging Spectrometer (HXIS), the X-ray 
Polychromator (XRP) and the Ultra- Violet Spectrometer 
and Polarimeter (UVSP). We concentrate on the obser- 



vations made with the XRP (Acton et al. 1980) and with 



the UVSP (Woodgate et al. 1980). In Peres et al. (1987) 



the light curves from the whole raster in the X-ray lines 
observed with the XRP were compared with the results of 
the previous version of the PH code. An extensive analy- 
sis and interpretation of the hot component of this flare 
as well as of the spatial heating distribution and evolu- 
tion was also included in that paper. During this flare 
the FCS instrument of the XRP registered spectroheli- 
ograms in six different resonance lines of some H-like ions 
or He-like: Oviii 18.97 A (T - 3 lO^K); Neix 13.45 A 
(T ~ 4 lO^K); Mgxi 9.17 A (T - 6.5 lO^K); Sixiii 6.65 
A (T - lO^K); Sxv 5.04 A (T ~ 1.55 lO^K); Fexxv 1.85 
A (T ~ 7 lO^K). These lines covered most of the coronal 
temperature range. 

There are no qualitative differences between the light 
curves synthesized with the new and the previous version 
of the PH code, except the disappearance of the numeri- 
cal noise which affected the results at the lowest coronal 
temperatures. Therefore we confirm that the previous sim- 
ulations carried out by Peres et al. (1987) are adequate to 
model the X-ray line light curves. 

Here we focus most of our attention on the UVSP 
data. This spectrometer had spectral resolution 1150-3600 
A and spatial resolution better than 2arcsec. The raster 
range was 256 x 256arcsec^. The UVSP collected images 
of this event with a field of view of 120 x 120arcsec^ and 
pixels of 10 X lOarcsec^ (i.e. ~ 7250 x 7250km2 on the 
Sun). Observations continued until the emitted flux in- 
creased considerably; since that moment (~17:00:06 UT), 
rasters (a complete raster lasted ^ 15s) were alternated to 



spectral scans in two different lines, which lasted 30s, with 
spectral resolution 0.3A. These two lines, simultaneously 
recorded in two separate detectors, are 

— Fexxi 1354.1 A line, formed at typical flaring plasma 
temperatures (T ^ 10 ''K); 

- Ov 1371.3 A fine, formed at TR temperatures (T 
2.2 lO^K). 

Both lines are partially blended with nearby lines 



( iFeldman et al. 1977|) . The Fexxi 1354.1 A hue partially 
blends with the chromospheric Ci 1355.844 A line; any- 
way this C I line is very narrow (~ 0.12 A), and, moreover, 
[Cheng et al. (1979) showed that, during the flare maxi- 



mum and decay phase, the intensity of the C i 1355.844 A 
line was ~ 20% of the Fexxi 1354.1 A line. The O v 1371.3 
A line instead appeared blended with a very narrow line 
at 1371.37 A (not yet identifled) during many flares ob- 
served by Skylab, which was not present in the quiescent 



phase (Feldman et al. 1977). 

The UVSP observed the rising phase in both lines. The 
flux in the O v 1371.3 A line, already high in the pre-flare 
phase, peaked at 17:02, i.e. at the same time as the peak 
observed in the hard X-rays. Then the emission in this 
line gradually decreased while the emission in the Fexxi 
1354.1 A line was still increasing (CP88). The light curve 
in the Ov 1371.3 A line was different from the other ones 
observed simultaneously, because it peaked twice in the 
pre-flare phase (at 16:58UT and at 16:59UT) and, more- 
over, the flux was higher where a fllament was observed 
in the Ha images. MacNeice et al. (1985)| wrote that ma- 
terial expulsions appeared in this fllament, at the velocity 
of 60km s~^ before the rising phase of the flare. This fll- 
ament disappeared after a short time and its disruption 
might be correlated to the increased emission in the Ov 
1371.3 A line. 

The light curves in the EUV lines recorded during this 
flare have been published in CP88. Cheng and Pallavicini 
distinguished three different zones in the flaring region, 
named as kernels Fl, F2 and LI. The two kernels Fl and 
F2 covered the two foot-points of the loop while kernel LI 
corresponded to the region between them, covering the top 
of the loop. The light curves in both lines were analyzed 
separately for each kernel as well as for the whole raster 
containing the entire flaring region, and are shown in Figs. 



From the model results we calculate the flux in a par- 
ticular line at a distance R from the Sun as 



F ^ 



2A 
47ri?2 



.G{T)ds, 



(1) 



where R is set equal to lAU in our calculations, A is the 
loop cross-section (assumed constant), G{T) is the emis- 
sion function tabulated by several authors for most ob- 
served lines, n is the hydrogen number density, Ue the 
electron number density and si and S2 are the extreme 
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loop points considered in the spatial integration; we mul- 
tiply the flux times 2 because the loop is supposed to be 
symmetric and therefore the integral is evaluated on half 
a loop length. 

As a first step the integral is calculated only in the 
loop volume corresponding to region Lf (two pixels). From 
Fig. 10 in CP88 we can guess that the loop was along the 
diagonal of kernel LI; assuming a semicircular geometry 
for the flaring coronal loop, and considering projection 
effects, the flux |l) is evaluated from si = 11 Mm to S2 — 
20 Mm along half a loop, i.e. over a distance 9 Mm from 
the loop top. 

Results for he Fe xxi 1354. 1 A are shown in Figs. ^ and 
^ and are obtained using the tables for the emission func- 
tion G{T) calculated by Monsignori Fossi (1994, private 
communication). The observed light curve appears simi- 
lar to those obtained for the X-ray lines with the FCS. 
This fact is not surprising as the FCS lines and the Fexxi 
1354.1 Aline originate in the same part of the atmosphere. 
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Fig. 2. Light curve observed in the Fexxi 1354.1 A line 
by the UVSP/SMM in kernel LI (sohd hue) from CP88. 
Also plotted are the light curves synthesized from hydro- 
dynamic models with the emissivity tables computed by 
Monsignori Fossi (1994, private communication). In order 
to allow a proper comparison, the emission is integrated 
around the top of the loop. 



The curve computed with simulation B (base heating) 
is closer to the observed one during the flare onset, and the 
shape of the light curve is well simulated for the first 400s. 
However the flux decreases too rapidly during the flare de- 
cay. On the other hand, with the heating at the top (simu- 
lation A), the curve is too low during the early phase, but 
it is high even well after the heating is switched off. The 
modeled light curve peaks later than the observed one, 
but follows the flare decay emission better, even though 
also model A predicts a more rapid decrease of the inten- 
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Fig. 3. Light curve in the Fexxi 1354.1 A line observed 
by the UVSP/SMM in the whole raster (sohd line) from 
CP88. Also plotted are the light curves synthesized with 
our code and with the emissivity tables computed by Mon- 
signori Fossi (1994, private communication) for both mod- 
els. The line emission is integrated along the whole loop. 



sity than observed, implying that further heating might 
be released during the flare decay in the EUV emitting 
plasma. 

In Fig. I the emission integrated over the whole loop 
is compared with the flux in the Fe xxi 1354.1 A line com- 
ing from the whole raster for both the simulations. The 
modeled curves are now very different from the observed 
one, as the total measured intensity is much higher than 
the calculated one in both cases. 

The observed emission in the Ov 1371.3 A line was 
brighter in both kernels Fl and F2, corresponding to the 
two foot-points of the loop, than elsewhere in the raster 
region. A chromosphcric filament was probably responsi- 
ble of the initially increased emission in kernel F2 during 
the pre-flare phase (see above). The observed flux in re- 
gion F2 is higher than in Fl and the two light curves are 
different (see CP88, Fig. 10), suggesting the presence of 
other emitting chromosphcric structures in region F2. For 
this reason we have compared our calculations only with 
the data of kernel Fl and not with those of F2. In kernel 
Fl there is high pre-flare emission, which could be related 
to the presence of other closed structures emitting in the 
Ov 1371.3 A line or in the blended unidentifled line men- 
tioned above. 

From our hydrodynamic calculations we simulate the 
light curve in this line using two different formulations for 
the emissivity G(T): 



1. the spectral emission function 



Landini & Monsignori Fossi (1990) 



tabulated by 
and following 



upgrades (thereafter LMF); 
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Fig. 4. The emission in the Ov 1371.3 A line of kernel 
Fl (covering one foot of the loop) is compared with the 
resuhs of our two numerical calculations and the LMF 
spectral synthesis model (top). In the plot at the bottom 
the same results are shown together with the flux emitted 
from the whole raster (CP88). 

2. the e missivity computed with Raymond's spectral 



code (Raymond fc Smith 1977 and following up- 



grades, henceforth RS). 

The two modeled light curves have the same profile but 
they differ by almost one order of magnitude. Discrepan- 
cies in the results obtained with different spec tral models 
have already pointed out by Pallavicini (1994) and consti- 
tute an important difficulty in comparing real data with 
the hydrodynamic calculations or any other modeling re- 
sults. Anyway, the discrepancy among spectral models 
does not alter significantly the shape of the light curve, 
but mostly the intensity. 

We have first compared the flux integrated over half 
a loop with the emission from kernel Fl; as a matter of 
fact, according to our model, only the TR contributes to 



Fig. 5. The emission of kernel Fl (covering one foot of the 
loop) is compared with the results of our two numerical 
calculations and the spectral synthesis model RS (top) . In 
the plot at the bottom the flux emitted from the whole 
raster (CP88) is compared with the same modeled curves. 



the emission in the Ov 1371.3 A line, and therefore, only 
the foot-points of the loop are bright in this wavelength. 
In Fig. ^ and in Fig. ^ we compare our numerical results 
(both simulations A and B) with the observed flux in ker- 
nel Fl (top panel). We also plot the light curves predicted 
by our hydrodynamic calculations and the measured flux 
coming from the whole raster (bottom panel). The raster 
flux is larger than that predicted with a single loop model, 
either if we assume the LMF spectral model or the RS 
model. This fact is not surprising because many other 
structures present in the raster might contribute to the 
emission, as for the Fexxi 1354.1 A line. Results obtained 
with the RS spectral model are in better agreement with 
the observations in kernel Fl. 

The Ov 1371.3 A line peaks earlier than the Fexxi 
1354.1 A line in our simulations as well as in the real ob- 
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servation but the O V 1371.3 A light curves of our models 
are qualitatively different from the data plotted in CP88, 
at least for the first 5 minutes. The model with the heating 
at the loop top and the RS spectral model approximate 
the decay light curve quite well. This could be said also 
for the bottom heated model for t < 700s. Then a thermal 
instability occurs (see long-dashed lines in Fig. |^ and in 
Fig. H). The rapid cooling of the atmosphere, mostly due 
to the increase of radiative losses as the temperature de- 
creases below lO^K, determines a significant peak of the 
emission which makes model B unlikely to describe this 
flare. 

5. Discussion and conclusions 

The X-ray emission of this flaring loop has been accu- 
rately modeled by Peres et al. (1987) using the hydrody- 
namic results of the PH code: in that work they obtained 
theoretical light curves very close to the observed ones. In 
this paper we compare the results of our numerical calcu- 
lations with both typical TR emission and coronal ones, 
simultaneously. 

From the analysis of all the modeled light curves in 
both the X-rays and EUV bands and the data registered 
during the event, we note that the top heated loop model 
gives results in better agreement with the coronal emission 
than the bottom heated loop model, during the entire flare 



ing coherently - dom inat ed the coronal em ission, as 



evolution, in agreement with Peres et al. (1987). On the 
basis of our new calculations we can add that this model 
also accounts for the observed Fexxi 1354.1 A line during 
a large fraction of the November 12, 1980 flare. We note 
that the computed light curves in the Fexxi 1354.1 A line 
have a shape similar to those of the X-ray lines formed 
in the same temperature range, and therefore originating 
from the same part of the atmosphere. 

Both our models predict Fexxi line flux close to the 
observed one at the time of the flare X-ray line peak, which 
occurred during the decay of the flare light curve in the ob- 
served O V 1371.3 A line. Also the top-heated model yields 
an already decaying OV light curve at the time of coronal 
lines flare maximum. However during the early flare the 
lower TR emission cannot be fitted with the same model 
which consistently yields lower emission than observed. 
This suggests that more plasma than modeled might have 
contributed to the Ov 1371.3 A line emission. 

In summary, the analysis of the hot coronal EUV 
emission (Fexxi 1354.1 A line) co nfirms the scenario 
derived with the previous study ( Peres et al. (1987) ) 
based on X-ray lines. The study of the Ov 1371.3 A 
line, instead, shows a rather different evolution and the 
relevant light curve is poorly fitted by the hydrody- 
namic loop model which instead works well for the hot 
coronal emission. We therefore confirm that the single 
loop model (albeit a schematic approximation of real- 
ity) works well for the corona of this fiare, and there- 
fore that a single loop - or a bundle of loops evolv- 



MacNeice et al. (1985)| and [Peres et al. (1987] suggested. 
The colder plasma emission, instead, appears dominated 
by other structures, probably some already present at an 
earlier phase, others evolving and reaching their maxi- 
mum emission before the Fexxi 1354.1 A line peak. Also 



Feldman (1983 )| , [Widing fc Cook (1987)| , [Poschek (1997) 



Mason (1998) observed that there is a mismatch between 
TR and coronal flare emission. 

It appears that a single loop model is not adequate 
to flt the lower TR lines simultaneously to the coronal 
lines (which are well fitted), at least for this flare. To this 
end more complex models may be required, such as multi- 
structure models. 
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